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In the target rest frame and at high energies, Drell-Yan (DY) dilepton production looks like 
bremsstrahlung of massive photons, rather than parton annihilation. The projectile quark is 
decomposed into a series of Fock states. Configurations with fixed transverse separations are 
interaction eigenstates for pp scattering. The DY cross section can then be expressed in terms 
of the same color dipole cross section as DIS. This approach is especially suitable to describe 
nuclear effects, since it allows to apply Glauber multiple scattering theory. We go beyond 
the Glauber eikonal approximation by taking into account transitions between interaction 
eigenstates. We calculate nuclear shadowing at large Feynman-x'F for DY in proton-nucleus 
collisions, compare to existing data from E772 and make predictions for RHIC. Nuclear effects 
on the transverse momentum distribution are also investigated. 



1 Introduction 

Although cross sections are Lorentz invariant, the partonic interpretation of the microscopic 
process depends on the reference frame. As pointed out by one of the authors, in the target 
rest frame DY dilepton production should be treated as bremsstrahlung, rather than parton 
annihilationS (see also 11). The space-time picture of the DY process in the target rest frame 
is illustrated in fig. |[ A quark (or an antiquark) from the projectile hadron radiates a virtual 
photon on impact on the target. The radiation can occur before or after the quark scatters off 
the target. Only the latter case is shown in fig. [I]. 

A salient feature of the rest frame picture of DY dilepton production is that at high energies 
and in impact parameter space the DY cross section can be formulated in terms of the same 
dipole cross section as low- a; Bj DIS. 

The color dipole approach to the DY process provides a convenient alternative to the well 
known parton model, in particular, it is especially appropriate to describe nuclear effectsu'0. 
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Figure 1: A quark (or an antiquark) 
inside the projectile hadron scatters 
off the target color field and radiates 
a massive photon, which subsequently 
decays into the lepton pair. The photon 
can also be radiated before the quark 
hits the target. Here, a is the longitu- 
dinal momentum fraction of the quark 
carried by the photon. 



2 DY dilepton production in pp scattering 

The cross section for radiation of a virtual jpJbxjton from a quark after scattering on a proton, 
can be written in factorized light-cone formH-HO, 

da(qp — > "f*X) 



, d p|* 7 i g (a,p)| a qq (x 2 ,ap), (1) 

similar to the case of DIS. Here, a qq is the cross section @ for scattering a gg-dipole off a proton 
which depends on the qq separation ap, where p is the photon-quark transverse separation and 
a is the fraction of the light-cone momentum of the initial quark taken away by the photon. We 
use the standard notation for the kinematical variables, x\ — x 2 = xf, t = M 2 / s = x\x 2 ^ where 
xf is the Feynman variable, s is the center of mass energy squared of the colliding protons and 
M is the dilepton mass. In (Q) T stands for transverse and L for longitudinal photons. 

The physical interpretation of (|]) is similar to the DIS case. The projectile quark is expanded 
in the interaction eigenstates. We keep here only the first eigenstate, 

\q) = V^\Qbare) + ^ L g \qi*) + ---, (2) 



where Z 2 is the wavefunction renormalization constant for fermions. In order to produce a new 
state the interaction must distinguish between the two Fock states, i.e. they have to interact 
differently. Since only the quarks interact in both Fock components the difference arises from 
their relative displacement in the transverse plane. If p is the transverse separation between the 
quark and the photon, the 7*g fluctuation has a center of gravity in the transverse plane which 
coincides with the impact parameter of the parent quark. The transverse separation between 
the photon and the center of gravity is (1 — a) p and the distance between the quark and the 
center of gravity is correspondingly ap. Therefore, the argument of a qq is ap. More discussion 
can be found ini. 

The transverse momentum distribution of DY pairs can also be expressed in terms of the 
dipole cross section!. The differential cross section is given by the Fourier integral 



da(qp -* 7*Y) 1 f j2 „ ^2 

din ad V (2vr) 2 



J d 2 p 1 d 2 p 2 exp[i<fj_ • (pi - p 2 )]^** q (a,p 1 )^ 1 * q (a,p2) 



x ^Wqq(x2,api) + o- qq (x 2 ,ap 2 ) - a qq (x 2 ,a(pi - p 2 ))} . (3) 

after integrating this expression over the transverse momentum q± of the photon, one obviously 
recovers (JlJ) . 

The LC wavefunctions can be calculated in perturbation theory and are well known II. The 
dipole cross section on the other hand is largely unknown. Only at small distances p it can 



be expressed in terms of the gluon density. However, several successful parameterizations exist 
in the literature, describing the entire function a qq -(x,p), without explicitly taking into account 
the QCD evolution of the gluon density. We use the parameterization by Golec-Biernat and 
Wiisthoffi for our calculations, fie. [2|. This parameterization vanishes oc p 2 at small distances, 
as implied by color transparency □ and levels off exponentially at large separations. The data 
in fig. |2] are quite well described without any iT-factor, which does not appear in this approach 
since higher order corrections are supposed to be parameterized in a qq -(x2,p)- 



10.0 



o 

.o 
c 

.9 6.0 



4.0 



S 2.0 



0.0 




7 8 
M/GeV 



10 




4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



x B =0.625 



7 8 
M/GeV 



10 



4.0 
3.5 
3.0 
2.5 - 
2.0 
1.5 
1.0 
0.5 
0.0 



x B =0.675 




7 8 
M/GeV 



10 



Figure 2: The points show the measured DY cross section in p H scattering from E772. The curves are calculated 
without any further fitting procedure. The solid curves are calculated at the same kinematics as the data points 
(center of mass energy ^/s — 38.8 GeV). The dashed curves are calculated for RHIC energy, yfs = 500 GeV. 



3 Proton-nucleus (pA) scattering 

Shadowing in DY is an interference phenomenon due to multiple scattering of the projectile quark 
inside the nucleus. In the target rest frame, where DY dilepton production is bremsstrahlung of 
massive photons, shadowing is the Landau-Pomeranchuk-Migdal (LPM) effect. These interfer- 
ences occur (fig. [l]), because photons radiated at different longitudinal coordinates z\ and z 2 are 
not independent of each other. Thus, the amplitudes have to be added coherently. Destructive 
interferences can occur only if the longitudinal distance z 2 — z\ is smaller than the so called 
coherence length l c , which is the time needed to distinguish between a quark and a quark with 
a 7* nearby. It is given by the uncertainty relation, 

I = _J_ = _J - ») M2 (4) 

c AP~ m N x 2 q\ + (1 - a)M 2 + a 2 m 2 



Here, AP~ is the light-cone energy denominator for the transition q — ► q^y* and q± is the relative 
transverse momentum of the 7*5 Fock state. For z\ — z 2 > l c , the radiations are independent of 
each other. 

An immediate consequence of this is that l c has to be larger than the mean distance between 
two scattering centers in the nucleus (~ 2 fm in the nuclear rest frame). Otherwise, the projectile 
quark could not scatter twice within the coherence length and no shadowing would be observed. 

We develop a Green function technique i, which allows one to resum all multiple scattering 
terms, similar to Glauber theory, and in addition treats the coherence length exactly. The for- 
malism is equivalent to the one proposed ina for the LPM effect in QED. Our general expression 
for the nuclear DY cross section reads 

dhia d\na 2 J J ^ J Zl J J 

- * 

* 7 * 3 (a, P2) PA (b, z 2 ) a qq - (x 2 , ap 2 ) G (p 2 , z 2 \pi,Zi) 
x p A (b,zi)a q g(x 2 ,api)^f 1 * q (a,p 1 ) . (5) 

The first term is just A times the single scattering cross section, where A is the nuclear 
mass number. The second term is the shadowing correction. The impact parameter is b and the 
nuclear density is pa- The Green function G describes, how the bremsstrahlungs-amplitude at 
zi interferes with the amplitude at z 2 . 

To make the meaning of Eq. I more clear, let us first consider a limiting case for G. In the 
simplest case, the coherence length, Eq. ||, is infinitely long and only the double scattering term 
is taken into account. Then G (p 2 , Z2 \ pi, z\) = 5^ 2 \pi — p 2 ) and one of the p integrations can 
be performed. The <5-function means that at very high energy (infinite coherence length) the 
transverse size of the 7*5 Fock-state does not vary during propagation through the nucleus, it is 
frozen due to Lorentz time dilatation. Furthermore, partonic configurations with fixed transverse 
separations in impact parameter space were identified a long time ago§ in QCD as interaction 
eigenstates. This is the reason, why we work in coordinate space. Namely, in coordinate space, 
all multiple scattering terms can be resummed and in the limit of infinite l c one obtains 

Gfrozen ^ ^ j ^ = j(2) _ ^ exp ^^Pl) ^ ^ (g) 

The frozen approximation is identical to eikonalization of the dipole cross section in Eq. (|l]). 
Thus, the impact parameter representation allows a very simple generalization from a proton to 
a nuclear target, provided the coherence length is infinitely long. 

At Fermilab fixed-target energies (y/s = 38.8 GeV for E772), this last condition is not 
fulfilled and one has to take a finite l c into account. The problem is however, that l c , Eq. |4], 
depends on the relative transverse momentum q±_ of the 7*g-fluctuation which is the conjugate 
variable to the size p of this Fock-state and therefore completely undefined in p-representation. 
The quantum mechanically correct way to treat the q\ in Eq. |I| is to represent it by a two- 
dimensional Laplacian Aj_ in p-space. The Green function which contains the correct, finite 
coherence length and resums all multiple scattering terms fulfills a two-dimensional Schrodinger 
equation with an imaginary potential, 

. d Aj_ (02) — ri 2 i . . s 

0Z2 Lhj q a (1 — a) 2 

= \5(z 2 -z 1 )5^(p2-pi). (7) 

For details of the derivation, we refer toi. 
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Figure 3: Comparison between calculations 
in the Green function technique and E772 
data at center of mass energy y/s = 38.8 
GeV. for shadowing in DY. The calcula- 
tions are performed at the mean values of 
the lepton pair mass. From the left to the 
right, these values are 5 GeV, 5.7 GeV, and 
6.5 GeV. 
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The imaginary potential accounts for all higher order scattering terms. The Laplacian implies 
that the Green function is no longer proportional to a ^-function. This means the size of the 7*g 
fluctuation is no longer constant during propagation through the nucleus. One can say that an 
eigenstate of size pi evolves to an eigenstate of size pi 7^ pi, so transitions between eigenstates 
occur. 

Calculations with Eqs. [5] and |7| are compared to E772 data§ in fig. §. Note that the 
coherence length l c at E772 energy becomes smaller than the nuclear radius. Shadowing vanishes 
as X2 approaches 0.1, because the coherence length becomes smaller than the mean internucleon 
separation. It is therefore important to have a correct description of a finite l c in this energy 
range. 

Nuclear effects on the (/^-differential cross section calculated at RHIC energy are shown in 
fig. ||. See@ for details of the calculation. The differential cross section is suppressed at small 
transverse momentum q± of the dilepton, where large values of p dominate. This suppression 
vanishes at intermediate q± ~ 2 GeV. In this region, one even observes an enhancement which 
reminds one of the Cronin effect. This enhancement is due to multiple scattering of the quark 
inside the nucleus. A nuclear target provides a larger momentum transfer than a proton target 
and harder fluctuations are freed, which leads to nuclear broadening. Note, that not the entire 
suppression at low q± is due to shadowing. Some of the dileptons missing at low q± reappear in 
this enhancement region. At very large transverse momentum nuclear effects vanish. 



4 Summary 

We express the DY cross section in terms of the cross section a q q for scattering a qq dipole off 
a proton. This is the same dipole cross section that appears in DIS. We can reasonably well 
describe low %2 DY data from pp collisions 13, without any free parameters and without a K 
factor. 

At very high energy, the dipole approach is easily extended to nuclear targets by eikonaliza- 
tion. At lower fixed target energies (E772) the frozen approximation is no longer valid, because 
the size of a Fock state varies during propagation through the nucleus. Therefore, transitions 
between interaction eigenstates (i.e. partonic configurations with fixed transverse separations) 
occur. 




Figure 4: Nuclear effects on the DY 
transverse momentum distribution at RHIC 
energy (y/s = 200 GeV) for carbon and lead. 
The figure shows the DY cross section ratio 
(da A /dx F dM 2 d 2 q ± )/(Ada p /dx F dM 2 d 2 q ± ) 
for dilepton mass M = 5 GeV and Feynman 
xf = 0.6. Here, q± — qr is the transverse 
momentum of the lepton pair. 
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We develop a Green function technique, which takes variations of the transverse size into 
account and resums all multiple scattering terms as well. Calculations with the Green function 
technique are in good agreement with DY shadowing data from E772. We have also calculated 
nuclear effects in the transverse momentum distribution of DY pairs at RHIC energy. The DY 
cross section is suppressed at low transverse momentum, but enhanced at intermediate q± ~ 2 
GeV. Nuclear effects vanish at very large q±. 
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